Asymptotics for a free-boundary model 
in price formation 



o . 

O . Maria del Mar Gonzalez * Maria Pia Gualdani 

^ I Univ. Politecnica de Catalunya UT Austin 

Q 



Abstract 

We study the asymptotics for large time of solutions to a one dimensional parabolic evolution 
equation with non-standard measure- valued right hand side, that involves derivatives of the solution 
computed at a free boundary point. The problem is a particular case of a mean-field free boundary 
' model proposed by Lasry-Lions on price formation and dynamic equilibria. 

, The main step in the proof is based on the fact that the free boundary disappears in the linearized 

problem, thus can be treated as a perturbation through semigroup theory. This requires a delicate 
choice for the function spaces since higher regularity is needed near the free boundary. We show 
, global existence for solutions with initial data in a small neighborhood of any equilibrium point, 

and exponential decay towards a stationary state. Moreover, the family of equilibria of the equation 
is stable, as follows from center manifold theory. 

_ Keywords: Reaction-diffusion, asymptotics, center manifold, free boundary, price formation 
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\0 ■ 1 Introduction 

We consider an idealized population of players consisting of two groups, namely one group of buyers 
of a certain good and one group of vendors of the same good. The two groups are described by two 
OO ! non-negative densities Jb, fv depending on {x,t) e M x R_|_. In the model, x denotes a possible value 

■ of the price and t the time. 

At a certain time t, the vendors would like to sell the good, and the function fv{x,t) describes the 
. ; ■ density of the vendors that are wiling to sell the good at price x. Meanwhile the buyers will try to get 

r> I the good at a cheaper price. The transaction takes place when the two groups agree on the price: we 

' denote by p{t) the agreement price. The price p(t) will be the highest price the buyers are willing to 

pay, and the lowest price the vendors agreed to sell the good. There exists a transaction cost, which is 
denoted by a positive constant a. When a buyer gets the good for the price p{t), the actual cost of his 
trade is p{t) + a, as well as the profit for the seller is p{t) — a. As a consequence, the buyer that got the 
good for the price p{t), will try in a later time to sell the good at least at the price p{t) + a and the 
vendor that sold the good for p{t) will try to get at a later time the same good for a price not higher 
than p{t) — a. Thus the parameter a introduces some friction in the system. 

The price p(t) results from a dynamical equilibrium between the two density functions. The ran- 
domness in the problem is measured by the diffusion coefficient of the two densities fv and fs, and is 
denoted by a parameter tr > 0. 

The above situation can be described by the following system of free-boundary evolution equations: 

/ %r-^& = Hi)^-=Pit)-a ifx<pit), t>0, 

\ fB{x,t)>Oifx <pit), /b(.t,<) =Offx>p(0, 
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together with 

^ - ^ 1^ = Ht)S.=pit)+a iix> p{t), t>Q, 2) 
fv{x,t)>0 ifx>p(<), fv{x,t) = Oiix<p{t), 

where 

^W = -T^(^W'^)-T^(^W,^). (1.3) 

The symbol S denotes the Dirac delta at the indicated point. The multiplier X{t) represents the number 
of transactions at time t, so (|1.3|) means that the flux of buyers which must be equal to the flux of 
vendors. The initial conditions 

/b(x,0) = fsix) and /v(.x,0) = (.x) 

are such that, for some pi in R, 

/^(x) > if X < p/, /b(x) = if X > 
(x) > if X > p/, /s(x) = if X <p/. 

The equation satisfy the property of conservation of mass. Indeed, both 

'P(t) i-+oo 

fv dx and / dx 

oc J p{t) 

remain constant for all < > 0. 



Equations (|l.ip - (|1.3p describe a mean-field model for the dynamical formation of the price of some 
good that has been very recently introduced in [13]. There the authors proved that, under suitable 
assumptions on smoothness and integrability on the initial data, there exists a unique smooth solution 
{fix,t),p{t)) for X e M of (fLT|) - (fL3)) . 

The most important question we are going to address here concerns the long time behavior of the 
system: will the good reach a stable price ( p{t) constant as < ^ oo ?) or will the price keep 
oscillating in time and never reach a stable value? 

We remark here that in a bounded interval with symmetric initial data then the solution remains 
symmetric for all times and asymptotics were proved by the authors in their previous work [lOj . However, 
the general case contains a new ingredient: a free boundary (see [1], for instance, for some background 
and examples on these type of problems). 

In this work we address the problem (|l.ip - (|1.3p in a bounded interval [—A,B], A, B > 0, for 
a < min{A/2, B/2}, with zero-Neumann boundary conditions. The aim is to show that if we start with 
an initial condition that is near a general equilibrium point in some suitable function space, then there 
exists a unique solution of (|l.ip - (|1.3p . that decays exponentially fast in time to a unique stationary 
state. In addition, in can be shown that the problem presents a two-dimensional family of equilibria, 
and that this family is stable. 

Although there is a well developed theory of semigroups and invariant manifolds for the study of 
evolution equations, the main novelty here is the fact that dynamical system arguments can be used for 
a problem that presents a free boundary. This is possible since we succeeded to treat the free boundary 
as a perturbation of the linearized problem. In fact, as we will see in the following sections, the free 
boundary disappears in the linearization and appears again in the nonlinear part of the problem as a 
term of lower order. 

This allows us to study the linearized operator with the classical semigroup theory and to get time 
estimates for the linear equation. Even though the linearized operator is a non standard one, we can 
explicitly compute its eigenvalues and corresponding eigenfunctions. Unfortunately the eigenfunctions 
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do not build an orthogonal basis with respect to the standard product in L^. This complicates the 
choice of functional spaces. 

Indeed, the choice of function spaces is a delicate step in the proof. They need to be big enough to 
allow delta functions in the equation, but on the other hand, higher regularity is needed near the free 
boundary. In order to give a more explicit characterization of those spaces, interpolation theory and 
pseudo-differential operators are needed. 

For simplicity of the notation, we rewrite the problem (|l.ip - ()1.3|) as the single equation 
df cr^ d'^f 



f{x,Q)^fi{x) ■xn[-A,B] 
/.(-AO)=/.(S,0) = 0, 



(1.4) 



where 



and 



f ■■= Jb- fv, fi = Ib- fv, and p(0) = pi. 

Moreover, for i > 0, 

f{p{t),t) = 0, /(.T,<) > for aU x € [~A,p{t)), f{x,t) < for all x e ip{t),B]. 

The unknowns are two quantities: f{x,t) (the solution) and p{t) (the free boundary). Assume, without 
loss of generality, that — 1. 



The problem satisfies important conservation laws. Indeed, condition (|1.3p implies preservation of 
goods and players: 

p(t) .B 

f dx = mi and — / f dx = m2, for all t > 0, 



'A Jp{t) 

with TOi :— J^'j^^fidx and m2 ■— — J^^ fi dx, mi,TO2 > 0. Note that mi + TO2 represents the total 
amount of players and m2 the amount of goods. 

There is an equivalent way to write the problem as a coupled system of equations. To see it, we 
differentiate the equation f{p{t),t) — 0, to obtain that 

Jx[P[t),t) 

Thus our problem is equivalent to the system 

ft = fxx - fx{p{t), t) [52;=p(t)_a - <^x=p(t)+a] ' 

fxx{p{t),t)) (1.6) 



p{t) - 



with boundary conditions 
and initial data 



Ix{p{t),t) ' 

/,(-A,t) = /,(S,t)=0, 
f[x,Q)^fi{x), p{Q)^pi. 
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To finish, let us remark that the present work (existence and decay results for small initial data) is 
the first step in the study of general asymptotics for this problem, which is an ongoing project of the 
authors with L. Chayes and I. Kim. In particular, in the recent paper [3], global existence in time and 
uniqueness for general initial data has been shown. The present paper is the only available information 
so far on the asymptotic behavior. 

Note that Chayes-Kim have recently studied an evolution Stefan problem in [4], [5] in relation to 
particle systems. Although apparently unrelated to our problem, it shares many of its features. 



2 Main results 

First we set up the problem in the interval [—^,5], for A,B > 0: we seek functions f{x,t) and 
p{t) e {-A + a,B - a), t e [0, +oo) that solve 

ft = fxx - fx{p{t),t) [Sx^p(^t)-a - ^x=p{t)+a] , 

f{p{t),t) = 0, f{x,t)>Oforxe[^A,p{t)), fix,t)<Oforx£{pit),B], (2.1) 
U-A,t)^ UB,t)^0, 

with initial data fi{x), pi satisfying 

h{pi) - 0, 



(2.2) 



fi{x) > for X e [-A,pi), fi{x) < for X e (p/,S], 
fj{~A) = fj{B) = 0, 
p{Q) ^pi. 

The initial datum has mass 

Pi i-B 

fidx = mi, — / fidx = m2, mi,m2>0. (2-3) 

-A J PI 

The equilibrium states of the above problem are well understood in Section [31 First, note that 
system (j2.ip - (l2.2p has infinitely many equilibrium points, i.e., functions which satisfy equation 

fxx - fx{P°) [Sx=p»-a - S,=pO+a] = 0, (2.4) 

where /"(p") — 0. More precisely 

r -X°{x~p"), (/-a,p" + a), 

f%x) = l -X°a, xe{p'' + a,B), (2.5) 

iA°a, xe{~A,p^^a), 

for A" > and / G {-A + a,B- a). 

We denote with foo the unique equilibrium point that is solution of the stationary problem (|2.ip - 
p.3p (which means the only equilibrium point /" that satisfies also the preservation of mass condition), 
and it is given by 

{-Aoo(x-poo), X £ {poo - a,poc + a), 
- Aoofl, X £ {poo +a,B), 

+ Xooa, X £ {-A,poo - a), 

where poo and Aoo are uniquely determined by the conservation of mass condition 

foe dx = 1711, - foodx^m2, 
-A Jp^ 
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and mi and m2 are the constants given by the initial datum in (j2.3p . 

We perform a perturbation argument. First of aU, we hnearize around any equihbrium point f'^ of 
(|2.4p . Let = 0. Assume, without loss of generality, that p° = 0. We seek solutions that are 

perturbations of this equilibrium, as 

where g is a solution of the problem 

gt = Lg + N{g). (2.6) 
Here L is the linearized operator (studied in Section 2]) given by 

Lg ■■= gxx - 9x{P°) [Sx=p°~a - Sx=pO+a] + 9{P°) [S'x^pO-a ~ ^x=p<'+a] > (2-7) 

and N{g) contains the rest of non-linear terms. Note that the free boundary p{t) does not appear in L. 
We first show that, apart from the zero eigenvalue, the rest of eigenvalues of L are strictly negative and 
isolated, with spectral gap given by 7 > 0. Thus the linear problem gt = Lg has a very nice solution 
coming from semigroup theory. 

In order to state our main theorem, we need to fix precisely the function spaces. The basic regularity 
we need to apply perturbation theory is L : X ^ X, with domain D{L) — Z, and N : Y X ior some 
spaces Z C Y C X . In particular, if we use fractional powers of operators, then we need Y — X" for 
some < a < 1. These fractional order spaces are simply the a-interpolation between X and Z . 

Naively, we would like to use the spaces 

Y ^ H'\-A,B), X ^ H''-^"{~A,B), Z ^ H''-^"+^{-A, B), 

for some < r < 1/2. We remark here that the functions we consider will have homogeneous Neumann 
boundary conditions, so they have a Fourier series expansion / ~ J2 fn^n- Thus Sobolev spaces 
H^{~A, B), r eM. can be defined through the norm 

ii/iiL- = E(i 

The space X is chosen so that L : X ^ X is well defined (note that the image of L contains derivatives 
of 5 functions). In fact for every < r < 1/2 one can find < a < 1 such that —3/2 > r — 2a. 

However, in order to get a Lipschitz estimate for the non-linear term of the form 

\\N{g)~N{g)\\x<C\\g-~g\\Y, 

we realize that higher regularity is needed near the free boundary p{t). More precisely we will need 
g € C^''^. Thus, we will ask the function to (locally) be in some for 3 < s < 7/2. Even in the 
symmetric case studied by the authors in [10^ . higher regularity at the free boundary was required. 

To handle this issue, we introduce a cutoff function around the point p'^ — 0. For the rest of 
the paper, we fix a small real number v > such that 4i/ << a. Assume that pj G (— t/, ;/). Let 
<j) G C^{—A, B) be a smooth cutoff function such that = in \—A, —2v\ U [2t^, B] and </> = 1 in [— i/, v]. 
Fix 0<r<l/2,3<s<7/2 and consider the space 

X -.^ {f cH'-'^°'{-A,B) : (l)f e H'-^''{-A,B)} , (2.8) 

with norm 
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On the other hand, we set 

Z {/ e iJ''"2"+2(-v4,B) : (f)f e H''-^°'+^{-A,B)^ , 

with norm 

ll/llz = + 

where (p is another cutoff function supported on (— 3i^, 3;^) such that = 1 on supp (f). Then we have 

Lemma 2.1. For any 3 < s < 7/2, there exist < r < 1/2 and < a < 1 satisfying all the above 
requirements such that L : Z ^ X is a bounded operator and Z is dense in X . 

The proof of Lemma lemma-bounded-operator is contained in [11) . that takes a close look at this 
kind of functional spaces. However, in order to make the present paper self contained, we rewrite the 
proof in the Appendix. 

For example, take s = 3.4, a — 0.9 and r — 0.2. Thus we can extend L to be an imbounded operator 
L : X ^ X with D{L) = Z. We should ask the initial data to be in the space 

Y '■= = {X, Z)a^2, 

given by the (real) interpolation between X and Z. Although we do not have a precise expression 
for Y, in the work [TT] we show a characterization of this interpolation space using pseudo-differential 
operators; in fact, it satisfies Y C Y C Y where 

Y:={feH^i-A,B) : cj,f e H%-A, B)} , 

and 

Y:={f€ H-{-A, B) : H%-A, B)} , 

for some smooth cutoff function (j) satisfying (f>= 1 on the support of cj). 

The main result of this work proves the existence of a unique solution f{x,t) for all time i > 0, if 
we start with any initial condition // close in norm of the function space 1^ to a general equilibrium 

state /°, say 

\\fl-f\\y<P, 

for some p > 0. We also establish exponential decay to a unique stationary state /oo (that might not 
be the same as /°, but it is uniquely determined from the masses mi, 7712). Of course, foo will be very 
near /°, but more significantly, the solution f{x,t) will not leave the neighborhood of size p, i.e., we 
have stability; center manifold theory provides a very elegant solution. 

Theorem 2.2. Fixed < 7 < 7, for any admissible equilibrium f^ defined in (12. 5p . there exist p > 0, 
C > such that if we start with initial data fi G Y satisfying 

\\fl-f^y<P, 

then there exists a unique solution f{x,t) of (|2.ip - (j2.3p for all time t > with /(x, 0) ~ fi{x) that 
satisfies 

/eCi([0,+cx3) 

and 

||/(-,i)-/oo||y <Ce-^*||//-/l^, (2.9) 

for all t > Q, where 7 is the spectral gap given in remark \4.3\ and foo is the unique stationary state that 
satisfies 

r f - r f r f - r f 

I foo — I fl, / Joo — / //■ 

J-A J — A Jpao Jpi 
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An additional interesting ingredient is to see whether the size p of the neighborhood can be taken 
uniformly on f'^: 

Theorem 2.3. The constant p in Theorem \2.2\ can he taken uniformly, independent of , as long as 
we restrict to the family of equilibria to the set A^, where 

{f'^ admissible : A*^ > x} ; 

for some x > 0. 

Some comments on the structure of the paper: Section [3] gives precise formulas for the family of 
equilibrium points. Section 3] computes the eigenvalues and eigenfunctions for the linearized operator 
L, while in Section [5] we apply semigroup theory to understand the linearized problem. The main 
steps in the proof of the theorem are contained in Section [SI first we give the necessary estimates on 
the non-linear part in order to have local existence; then we show existence and stability of a center 
manifold so that we have global existence and decay for small solutions, and finally we comment on the 
stability. 



3 The stationary state 

First we review the construction of the stationary solution from [13'. Note that, each equilibrium point 
is determined by two quantities (the masses mi and m2). 



Lemma 3.1. Given any two positive real numbers mi,m2, the system Il2.1\) - i[273\) has a unique station- 
ary solution fac (x) , Poo */ only if 



2A + 2B-ia ~ m2 
Proof: Any equilibrium point of (|2.ip - (l2.2p must satisfy 

fxx [5x=p«-a - 5x=pO+c] m B] X R+ 

/,"(-A) = /°(B)=0 

/" > if .T e /° < if X G {p°,B). 

The solution is a piecewise linear function given by 

{0 if .TG (0,39° -a) 
/°|pO ifa:e (p"-a,p° + a) 
if a; G (/ + a,yl). 

The unique stationary state {foe, Poo) is computed with the above formula, if we impose the conservation 
of mass property, 

/Poo i-B 
foo dx, 7712 := - / foodx. 
-A Jp^ 

Since /oo(Poo) = and Aoc = ~dxfoo{Poo) > 0, then 

TOl = Aoofl {Poo - f + ^) I rn2 = AooO (i? - Poo - f ) , 

thus, 

mi _ Poo - ^ + A 
m2 -B - Poo - f ' 
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Note that the quotient nii jm-i is an increasing function of , as expected. From here we get 

_ -a{mi - 1112) - 2Am2 + 2Bmi _ mi + m2 , , 

Pec — ■ T , '^oo — —7 ; ; j^- (."J-^j 

2(mi+m2) a[—a + A + B) 

To finish, the condition that Poo G ^ + a, S — a) is equivalent to (|3.ip . □ 

Remark: We we say that an equihbrium Z*', f^{p^) = as above is admissible if it satisfies (|3.ip . Note 
that the set of admissible equilibria is a smooth family parameterized by mi,TO2, or by p'^,\^, and it 
has dimension two. 



4 Linear stability 

In this section we look at the eigenvalues of the problem in the interval \—A,B]. Let us assume that 
we are in a situation where an admissible stationary state /° exists and satisfies p° — 0. Remark that 
if m-i > 77i2, then A > B, and analogously, if mi < TO2, then B < A. On the other hand, toi — m2, 
then A = B, and the steady state is an odd function, defined on the interval [—A, A]. 

First, we linearize the equation around the stationary state f°,p^. If / is a small perturbation of 
f^, then we can see that it must have a unique root near p", call it p{t). This perturbation must satisfy 

ft = fxx - fx{p(t),t) [(5a;=p(t)_Q - 5x=p(t)+a\ 

< f{pit),t)^0 
J,{-A,t) = UB,t)=0. 

We write f = f'^ + eg, p = p'^ + eq, g ~ g{x, t), q — q{t), differentiate in e and set e = 0. We obtain 

gt = gxx - {fLiP%it) + 9x{p",t)} [<5,=po„, - <5,^po+„] - [S',^po_, ~ <5Upo+J q{t), 

together with the boundary conditions gx{—A,t) ~ gx{B,t) — 0. Use that = 0, and set 

A" = -fSip°) so that 



i 9t= gxx - gx{p°, t) [4=pO-a - 4=pO+a] + A° [^UpO_„ - (5^=pO+„] 

\g,i-A,t)=g,{B,t) = 0. 
On the other hand, we linearize the condition f{p{t), t) = and this gives 

so when we substitute the above formula in (|4.ip , we get the linear equation 

J St = gxx - .9x(p°, t) [S^=pa_a - 5^=pO+a\ + 9iP°) [S'x=pO-a " ^'x=pO+a\ > 

\g,{^A,t)^g,{B,t)^0. 
Note that Lemma [O] assures that the operator L is well defined. 



(4.1) 



(4.2) 



In the following, we show that the pointwise spectrum of the operator L consists only of real, 
non-positive eigenvalues. 

Proposition 4.1. Consider the operator 

Lg g,, - g,(0) [,5,=_, - 4=a] + 5(0) [<5U-a " 4=a] , 

defined on the space X with x £ [— A, _B]. Its pointwise spectrum can be described by 
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• Zero is an eigenvalue with eigenspace of dimension two. Two linearly independent eigenfunctions 
are given by 

X, X E {~a,a) ( 1, x £ {—aja) 

9o{x) = { a, xe{a,B), ho{x) = < 2, xe{a,B) 



xe{~A,~'a) I 2, xe{~A,-a). 



• The rest of the eigenvalues must be of the form 

2 . m: nn rnr 

pi = -a , for a= — ,— ,— , for some n e Z. 

a 2A — a 2B — a 

Every eigenspace has finite dimension. 

Proof: Ideas come from the previous work of the authors [TUj on the symmetric case. Let /i € C, and 
g & X such that Lg = /ig, with boundary conditions 

g^{-A)=g^{B)=0. (4.3) 

More precisely, 

gxx - 9x{0) [5x=-a - <5:c=a] + g(0) - 5'^^^] = fig. 

We see that the function g must be discontinuous at —a and +a with jump, 



5(a+)-g(a-) = .g(0), 
g{-a+)~gi-a-) = ~giO), 

and that its derivative too, with jump 

g'{a+)-g'{a-)^-g'{0), 
g'{^a+)-g\-a-)^g'{0). 



(4.4) 



(4.5) 



It is clear that, except at the points x — ±a, the function g must satisfy gxx — fig. Thus, our main 
trick is to find our eigenfunctions in a similar way as one finds eigenvalues and eigenfunctions for the 
Laplacian using Fourier series. 

Let a = e C, a 7^ 0. We set 

5i(a;) = cie"^ + C2e""^ in (-a, a), 

52(x) = die"" + d2e-"" in(a,S), 
53(2:) = eie"^ + 626""^ in(-A, -a). 



This g must still satisfy the boundary conditions (|4.3p and the matching conditions (|4.4p - (|4.5p as given 
above. When we impose (|4.3p . it is easy to see that g can be written simply as 



gi(a;) = cie°^ + C2e in (-a, a), 

<?2(a;) = rf[e-2"S+ax^g-a.j m{a,B), 

^3(2;) = e[e2"'4+"^ + e-"^] in(-A,-a). 
The matching conditions at a; = a imply that 

J d [e'^»B+ca _^ g-aaj _^ ^^^aa _ 1] + ^2 [-e""" - l] = 0, 

\ d [e-2"-B+"° _ e-"''] + ci [-e"° + 1] + C2 [e""" - l] = 0. 



9 



(4.6) 



When we add and substract the two equations above we get: 

cie"" + C2 + de-2"S+aa ^ 

ci + C2e-"° - de-"" = 0. 

On the other hand, we can substitute the matching conditions at x = —a to obtain: 

ci [e-"° + 1] + C2 [e""^ + 1] - e [e2"^-"" + e"*^] = 0, 
ci [e-"° - l] + C2 [-e"" + 1] - e [e^"^-"" - e"°] = 0. 

Then, adding and substracting the previous two equations: 

ci + C2e"" - ee"° = 0. ^ ' 

Then we see that (|4.6p together with (|4.7p become an homogeneous hnear system of four equations 
with four unknowns. Performing row reduction in the system, we can compute the determinant of the 
coefRcient matrix, and indeed, it is a (non-zero) multiple of 

It is clear that system (|4.6p - (|4.7p only has the trivial solution unless one of the three factors above 
vanishes. In this case, we must have that the real part of a is zero, and the imaginary part must take 
the values: 

X a) = — ,— ,— , 71 eZ. 

a ZA — a ZB — a 

The eigenvalues must necessarily be of the form ji — , which are all real and negative. The eigenspace 
for each eigenvalue is finite dimensional and depends on the number of solutions of the linear system 
above. 

To complete the proof of the proposition, we need to check the zero eigenspace. Then we seek 
functions g such that 

g" - g'{0) [S.,= -a - S^=a] + 9{Q) K.-. - 4=a] - 0. 

In particular only the following functions go and ho (and a linear combination of those) will be suitable 
solutions: 

{.go(0) = and (14. 5|) is satisfied, 
h'^{Q) = and is satisfied. 

The corresponding zero eigenspace is given by a linear combination of go and ho'. 

X e (—a. a) r 1, X {—a, a) 

go{x) = { a, X £ (a.B) , ho{x) = < 2, x e{a,B) . 

a;e(-A, -a) [2, x£{-A,-a) 

It is an easy computation to check that go and ho are linearly independent. 

□ 

It is clear from the method above how to compute the eigenfunctions. For completeness here, we 
show the explicit calculations in the case A = B = 1. 

Proposition 4.2. Consider the operator 

Lg := gxx - gx{0) [5x = -a - 5x=a] + g(0) [5'x = -a - ^'x=a] ) 

defined on the space X with x G [—1, 1]. All the non-zero eigenvalues of L are given by: 
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1. For all n such that 



we have 



2n 1 
T~ 2 ^ 



-{an) for Un := — 



with corresponding one dimensional eigenspace with eigenfunction given by 

sin(aa;) x G {—a, a) 

9n{x) = "( X e (a, 1) 

0, xe {-I, -a). 



2. For all n such that 



we have 



(2n + 1) 1 
a 2 



IJ-n = -{a„f for a„ := 



(2n + l)7r 



with corresponding two dimensional eigenspace given by the linear combination of the following 
functions 



gn{x) = < 



sm{ax) 
2 sva.{ax) 



X e (—a, a) 
X e (a, 1) 



2sva.{ax), a;e(— 1,— a). 



and 



hn{x) = < 



cos(Q;a;) 



0, 



X e (—a, a) 
X G (a, 1) 
X e (—1, —a). 



3. For all n such that 



we have 



(2n + l) 1 

a 2^' 



IJ-n = -{a„f for a„ := 



(2n + l)7r 



with corresponding one dimensional eigenspace generated by 

gn{x) = < 



cos(Q:a;) X € (—a, a) 

X e {a, 1) 

_0, xG{-l,-a). 



4- For all n such that 



have 



2n 1 ^ „ , 2n ^ „ 

^ Z, and — ^ Z, 

a 2 a 



-(q:„)^ /or a„ := 



2mT 



with corresponding one dimensional eigenspace generated by 
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5. For all n such that 



we have 



sin(a3;) x G {—a, a) 

9n{x) = { X e (a, 1) 

0, xG{-l,-a). 



S: ana — e 

a 2 a 



IJ'n = -{anf for an := 



2mr 



with corresponding two dimensional eigenspace generated by 



and 



sin(ax) 
gn{x) = { 
0, 



X e {—a, a) 
X e (a, 1) 
X e (-1,-a), 



cos(aa;) 
hn{x) = i 2cos(aa;) 



X € (—a, a) 
X G (a, 1) 
2cos(aa;), x&{—l,—a). 



6. For all n such that 



we have 



2n 1 ^ , 2n ^ ^ 
— - ^ Z, and a f. Z, 



2- a 2 



2-a 



^„ = -{anf for an 



2mr 
2-a 



with corresponding two dimensional eigenspace generated by 



s'm(a{l — a)) . , , cosfafl — a)) , . 

sm(aa;) + — -^-^ — -—^cos{ax), 



gn{x) 



2(1 — cos(Q;a)) 
cos(q;(1 — a;)), 
0, 



2(1 + cos(aa)) 



and 



sin(a(l — a)) , , cos(q;(1 — aX) 

^ ^ " s\a{ax) + — — ^cos(aa;), 



hn{x) = < 



2(1 + cos(aa)) 



2(1 — cos(aa)) 
0, 

cos(q!(1 + x)), 

Moreover, for all n> 1, the eigenf unctions have zero mass, i.e., 



x e (—a, a) 

X e (a, 1) 
X e (-1,-a), 



X € (—a, a) 

X S (a, 1) 
X G (—1, —a). 



/o M !■ 

gn dx = / hn dx = 0, / gn dx = j hn dx = 
-1 J-i Jo 



(4.8) 
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Proof: We have seen that all the eigenvalues are of the form /i = — < 0. The construction of the 
eigenfunction g is done piecewise in each of these three intervals: (—a, a), (— 1,— a), and (a, 1), using 
Fourier series. We seek 

gi{x) ~ ci sin(Q!a;) + C2 cos(Q!a;) in (—a, a), 

g2{x) — d\ sin(aa;) + d2 cos(aa;) in (a, 1), 

33(2;) — ei sin(aa;) + 62 cos(aa;) in (—1, —a), 



with boundary conditions 

.93(-l) = .92(l) = 0. 
The matching conditions (|4.4p . (|4.5p are rewritten as 

ff2(a) -51(a) =51(0), 
ffi(-a) - .93(-a) = -.91(0), 
52(a)-5;(a) = -g'i(0), 

ffi(-a) -53(-«) = 3i(0)- 
The zero Neumann boundary conditions (|4.9p give that 

di cos(a) — (i2 sin(a) = 0, 
ei cos(a) + 62 sin(a) = 0. 

Consider first the case cos(a) = 0. 

This implies ^2 = 62 = and from the matching conditions (|4.10p we get the system 

d\ sin(aa) — c\ sin(aa) — ci cos(aa) = C2, 
d\ cos(aa) — c\ cos(aa) + C2 sin(aa) = — ci, 
-ci sin(Q!a) + C2 cos(aa) + e\ sin(Q!a) = — C2, 
Ci cos(Q!a) + C2 sin(Q!a) — e\ cos(aa) = c\. 



(4.9) 



(4.10) 



(4.11) 



Consider now the case sin(aa) 7^ 0: eliminating the constant c\ from the first and second equation 
in (|4.1ip . we get 

d\ sin(Q!a) = 0, 

which is satisfied only if d\ — 0. With a similar computations in the third and fourth equations, we get 
ei sin(aa) = which is satisfied only for e\ = 0. Therefore system (|4.11l) reduces to 

— ci sin(aa) — C2 cos(aa) = C2, 
-ci cos(aa) + C2 sin(aa) = — ci, 
-ci sin(aa) + C2 cos(aa) = — C2, 
ci cos(aa) + C2 sin(Q!a) = ci, 

which is the unique solution ci = C2 = 0. 

Therefore consider now a such that sin(Q!a) = 0: this yields to 

— C2 cos(aa) = C2, 
di cos(aa) — ci cos(aa) — — ci, 
C2 cos (aa) = — C2, 
ci cos(aa) — ei cos(aa) = ci, 



which has solutions: 



di = ei = C2 = 0, for cos(aa) — 1, 
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and 

di = 2ci, ei = 2ci, for cos{aa) = — 1. 

This proves Assertion (1) and (2). 
Consider now the case cos(a) ^ 0. 
We can write 

sin(a) sin(a) 

di=d2 ei = -62 

cos(a) cos(a) 

which imphes that from the matching conditions (|4.10[) we get the following system 

d2 cos(q!(1 — a)) — ci sin(Q!a) — C2 cos(aa) — C2, 
d2 sin(a(l — a)) — ci cos(Q!a) + C2 sin(Q!a) = — ci, 
-Ci sin(Q;a) + C2 cos(aa) — 62 cos(a(l — a)) ~ — C2, 
Ci cos(aa) + C2 sin(aa) + 62 sin(Q;(l — a)) — ci. 



(4.12) 



Suppose sin(Q;a) = (we remind the reader that we are under the assumption cos(q;) 7^ 0): the 
above system reduces to 



(4.13) 



d2 cos(Q!a) cos(q;) — C2 cos(Q!a) = C2, 
c?2 sin(Q;) cos(Q!a) — ci cos(aa) = — ci, 

C2 C0s(Q;a) — 62 COs(q;) COs(Q!a) = — C2, 

ci cos(aa) + 62 sin(a) cos(Q;a) = ci. 
If cos(Q;a) = — 1, system (|4.13|) has a unique solution 

^2 = 62 = Ci = 0, 

and Assertion (3) is proven. 

Otherwise, if cos(aa) = 1, we separate the case when sin(a) = and sin(a) ^ 0. 
For sin(a) 7^ system (j4.13p has solution d2 = 62 ^ C2 — 0, which proves Assertion (4). 
For sin(Q?) — system (|4.13l) has a unique solution 

d2 = 2c2,e2 = 2c2, 

and Assertion (5) follows. 

Suppose sin(aa) 7^ and again cos(q!) 7^ 0. Eliminating ci from the first and second equation of 
(|4T^ we get 

^2 cos(a(l — a)) = ^2 cos(a). 

On the other hand, eliminating C2 gives d2Sm(a) = —d2sm(a{l — a)), which implies, that d2 ~ 
or a — l^^. Similar computations using the third and fourth equation of (|4.12p yield to 62 = or 
a = However, we can check that c?2 = 62 = does not produce any non-trivial solution. For 

a = combining the first and third equations of (|4.12p we get that 



d2 cos(a(l — a)) + 62 cos(a(l — a)) ~ 2c2(l + cos(Q;a)), 

which implies 

1 / ^^ ^2+62 

62 = - cos(q;(1 — a))- 



2 "(l + cos(aa)) 

Similarly, from second and fourth equation of (|4.12D we get 

61 = -sm(6i;(l - a))- 



2 ^ ^ ^^(l + cos(Q;a))' 
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Assertion (6) is proven. 

Summarizing, the structure of the eigenvalues foUows the foUowing scheme: 

sin(Q;a) 7^ no solutions 
cos(Q;a) = 1 



cos a 



cos(a) ^ < 



sin(aa) — ■ 



sin(Q;a) — 



cos(Q!a) — 1 < 
sin(aa) 7^ (Assertion 6). 



(Assertion 1) 
cos(Q;a) — 1 (Assertion 2), 

cos(aa) = — 1 (Assertion 3) 

sin(a) = (Assertion 5) 
sin(Q;) 7^ (Assertion 4), 



Remark 4.3. We write the spectral gap as 




5 Semigroup theory 

In the following, we study the properties of the linear part of our equation (|2.6p . given by gt = Lg, 
where the operator L is defined as in (j2.7p . In particular we will show that it generates an analytic 
semigroup, together with time dependent decay estimates. We close the section with a very explicit 
characterization of the ker (i), where 

ker(L) = span(go,/io), 
with go and ho defined as in Proposition [521 

To begin with, we review some important concepts from semigroup theory. Standard references are 
the first chapter of the book [12] and also [M]. Let X be a Banach space with norm \\-\\x and let L be 
a linear operator on X that has domain and range in X. Denote by cr{L) C C its spectrum. 

We say that L is a sectorial operator if it is a closed densely defined operator such that, for some 
(p £ (0,7r/2) and some M > 1 and real a, the sector 

Sa,4, := {A e C : < |arg(A- a)| <n,X^a} 

is in the resolvent set of L and 

\\{XI - L)-^\\^ <M/\\-a\ for aU A e S'a,0. 

An analytic semigroup on a Banach space A is a family of continuous linear operators on A, {T{t)}t>o, 
satisfying 

1. r(0) = /, T{t)T{s) = T{t + s) for i > 0, s > 0. 

2. T{t)x -^xast^O+, for each x € X. 

3. i — > T{t)x is real analytic on < i < 00 for each x £ X. 
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The infinitesimal generator L of this semigroup is defined by 



+0+ 



its domain D{L) consisting of all x G X for which this limit exists. And viceversa, it is well known 
that (see Theorem 1.3.4 in fT2], for instance) if L is a sectorial operator, then —L is the infinitesimal 
generator of an analytic semigroup {T{t) = e^*-^}t>o; this semigroup gives the solution of the ODE 
gt + Lg = 0. Moreover, if dia{L) > 7, then for t > 0, we have the bounds 



|e-^*||^ <Ce-^*, ||ie-^*||^ < ^e"^*, 



for some constant C. This implies in particular that Range (e ) C D{L). 

Suppose that L is a sectorial operator and ^a{L) > 0. Then for any a > we define the (— a)- 
fractional power as 

1 f°° 
L-" = ^-r / t^'-^e-^'dt. 



r(a) 

The operator L~" is a bounded linear operator on X which is one-to-one. For a > 0, L" is defined to 
be the inverse of L~", D{L°') = R{L^°'). L" is a closed, densely defined operator. We define 

X" := D{L"), 

with the graph norm 

It is seen in the last section that the space X" is precisely the a-interpolation between the spaces X 
and D{L). 

If the spectrum of the operator L has a good structure, then we can split the space X. A set 
(T C cr{L) U {00} —: (y{L) is a spectral set if both a and a{L)\a are closed in the extended plane 
CU {00}. 

Theorem 5.1 (Theorem 1.5.2. in |12|). Suppose L is a closed linear operator in X and suppose ci 
is a bounded spectral set, and a2 ~ a{L)\ai so 02 U {00} is another spectral set. Let Ei, E2 he the 
projections associated with these spectral sets, and Xj = Ej{X), j=l,2. Then X — Xi © X2, the Xj 
are invariant under L, and if Lj is the restriction of L to Xj, then 

Li : Xi Xi is bounded, cr{Li) ~ (Ti, 

D{L2) = D{L) n X2 and (7(^2) = <J2- 

The following lemmas will allow us to compare our operator L to the standard Laplacian. The first 
one will be used in order to prove that our linear operator L is sectorial, while the second gives that 
the spectrum of L := A -I- i? consists only of eigenvalues. 

Lemma 5.2 (Theorem 1.4.5 in |12j). If M is a sectorial operator with ^a{M) > and if B is a linear 
operator such that BM^°' is bounded on X for some < a < 1, then M + B is sectorial. 

Lemma 5.3 (Chapter II, section 5.13, and Chapter V, Corollary 1.15 in 8 ). Let M be an operator 
having compact resolvent and let B be a bounded operator on X such that the resolvent set of M + B 
is non-empty. Then M + B has compact resolvent. As a consequence, the spectrum of M -\- B consists 
only of eigenvalues. 
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The next lemma is used in order to define the domain of a fractional power of an operator X" {— Y) , 
and to get a time decay estimate for the semigroup. This space is important because it is the space 
of admissible initial conditions; the lemma tells us that it is enough to consider fractional powers of 
the Laplacian to define it. However, a more explicit characterization of Y in terms of usual norms is 
desirable - this is done in [11] . 

Lemma 5.4 (Theorem 1.5.4. in [12]'). Suppose L is a sectorial operator, ai a bounded spectral set for 
L, (72 = a{L)\ai, 3?(T2 > 7 and X = Xi ® X2 is the corresponding decomposition. Assume also that M 
is a sectorial operator with D{M) = D{L), 5Rct(M) > 0, (M — L)M~" is bounded for some < a < 1. 
Then using the norm \\x\\^a. := ||M"a;||, < a < 1, /or x G n D{M") and t > 0, 

for some positive constant Ci . 

Consider now our linear problem 

9t = Lg, 

g{x,0)^gi, (5.1) 
[g,i-A,t)=g,iB,t) = 0. 

where L : X X is defined as in (|2.7p . and the initial condition gi G Y. Using the results recalled 
above, we are able to prove now the next proposition: 

Proposition 5.5. Let Ui ~ {0} and a-i be the rest of non-zero eigenvalues of L, defined as in (|2.7p . 

Let El, E2 be the projections associated with these spectral sets, and Xj = Ej(X), j—1,2. Then we 
have the splitting X = Xi ® X2, where the Xj are invariant under L, and if Lj is the restriction of L 
to Xj, then 

Li : Xi Xi is bounded , (t(Li) = tJi, 
D{L2) = D{L) n X2 and (7(2.2) = (72- 

In our case, Ei — ker(Z/), Li = 0. Moreover, —L2 is sectorial generates an analytic semigroup T{t) = 
e^^*. This semigroup satisfies the estimates 



\\Tit)x\\y<Ci\\x\\^t-<^e 
\\T{t)x\\y < CiMye-^' 



(5.2) 



for any < 7 < 7. 



Proof: Proposition 14.21 computes the spectrum of L. Apart from the zero eigenvalue (with eigenspace 
of dimension two), the operator —L has a discrete and positive spectrum. Thus Theorem 15.11 splits the 
space X: let CTi = {0} and a2 be the rest of eigenvalues. Then X = Xi (B X2 where Xi = ker(L), that 
satisfies dim(Xi) = 2. 

Next, we write — L2 = —A + B\^^ where 

B = (.9x(p°) [Sx=pO-a - 4=pO+a] " 9iP°) [S'x^pO-a " ^x=pO+a\) ■ 



However, — A|^^ contains only positive eigenvalues. Thus Lemma 15.21 applied to AI = — and B as 

above gives that —L2 is sectorial, once we note that _D(— A? ) C D{B). In fact, since dI—A? ) C 

1x2 

D(B), closed graph theorem implies that BA7°' is bounded. 

1x2 
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Next, we use Lemma [5.41 for M ~ ^^ix, ' ^'^'^ /3 = a G (0, 1), so 
for some Ci positive constant. 



□ 



The solution of equation (|5.ip can be written now as 

g[x, t) = cigo + diho + e^^*-{gi), 
with cigo+djho G Xi for some c/, c?/ g R that depend only on the initial condition gi, and e^^*^{gi) G X2. 

Although it is not usually known, in our situation we are able to give a more explicit characterization 
of Xi — keri. It will be needed later. 

Lemma 5.6. Given any g G X, g — gi + g2, gi S Xi, i = 1,2, then in the basis of eigenf unctions 
{50, /lo} given in Proposition \4.S\ we can write 

91 = ego + dho G Xi, 

for 

1 



( 94V ^r,Ai-a + 2A)h[g]-{-a + 2B)h[g]}, 
a[a — ZA)[a — IB) 

d = ^ {(-4 + aB)h[g] - (4 - aA)h[g]] 



(5.3) 



a{a ~ 2A){a - 2B) 
where Ii [g] and I2 [g] are defined as 



g dx, hig] := / g dx. 

-A Jo 



Proof: We know that g = ego + dho + <72, for 52 G X2 for some c, d. If we integrate this expression we 
obtain 

B pB pB i-B 

g dx — c I go dx + d ho dx + g2 dx, 
"'0 "'0 "'0 

pO M M 

g dx = c I go dx + d ho dx + g2 dx. 

-A J-A J-A J-A 

However, because g2 G X2 (see (|4.8p ) we know that 

pB 

g2 dx = / (72 dx = 0, 
'-A Jq 

thus we conclude from the previous calculation that c, d must solve 



l2[g] = c [ 
Jo 

li[g]^c [ 
J —J 



B pB 

go dx + d ho dx, 
Jo 
no 

go dx + d ho dx. 



-A J-A 
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The integrals of go, ho can be explicitly computed; we obtain a system of two equations and two 
unknowns c, d 

l2[g] = c (-^ + aS) + d{2B - a), 

Ii[g] = c - aA) + d{-a + 2 A), 

whose solution is straightforward. □ 

Remark: Note that when wc apply the previous result to our case, the functions Ii[g] and l2[g] will 
be time depending functions, more precisely 

/PI l-B rp(t) 

fidx- fdx- [f + g] dx, 
-A Ja Jo 

l-B fB ML) 

h[g] = fidx- fdx+ (f + g) dx, 

Jpr Jo Jo 

where pi is the unique zero of the initial data // and p(t) the unique zero of the function f{-,t) — 



6 Non-linear theory 

Let be an equihbrium point as specified by the hypothesis in the theorem. We try to construct 

a solution of equation 

ft = fxx - fx{p{t)) [^p{t)-a - 5p{t)+a\ 

as a perturbation of this stationary state, i.e., 

f{x,t) = f{x)+g{x,t). 

Consider the space V of functions / : [—A, S] ^ M that have only one root p G (p° — + v) and 
that are positive in [—A,p), negative in {p,B\. We will allow perturbations g such that g{-,t) for each 
fixed time belongs to the space 

Up = {geY : \\g\\y<p, f+geV], (6.1) 
for some p. It is essential that this neighborhood is open and uniform in some sense: 
Lemma 6.1. Given p > 0, and < iv < p, if u) < min |A°j^, A"}, the neighborhood 

{geY : \\g\\y<uj} 

is contained in Up. 

Proof: Assume that ||.g||y < w as specified in the hypothesis of the lemma. Let us prove that / = f^+g 
has an unique zero in the interval {—u,!/). First of ah note that / is strictly monotone decreasing in 
the interval {—y,v) since 

f,{x) = f^{x)+g,{x) = -X''+gxix) < -A° +a; < 0, 

for LO < X°. Moreover 

/(-z/) = |A°|i/ + g{-u) > \X°\iy - w > 0, f{v) = -\X°\iy + g{v) < -\X°\v + w < 0. 
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This implies that / has an unique zero in the interval {—v, v). Let us call this unique zero p. It holds 
that —fxip) > X''^ — w > 0. A similar argument gives that / is positive in [— U {p, B]. □ 

Now we relate the nonlinear problem to the linear one, in order to obtain a formula for the non- 
linear part. Let f{x,t), p{t) be a solution to our problem. We write f{x,t) — f^{x) + g{x,t), and 
p{t) = p° + q(t). Assume that q G (— j/, i'). As in the calculation in Section [H we expand for each time 
t. Thus 

= t) = + + t) + (6.2) 

where 

Riit)=g{pit),t)-gip'',t). 

But — and = — A°. For simplicity in the notation, we drop the dependence in t. Then 

from (|6.2p we can write 

9= Yo • (6-3) 

On the other hand, we expand 

Upit),t) = + ff,(/,<) + i?2 = - A" + g,(/,<) + i?2, (6.4) 

where 
and also 

When we substitute (|6.3p . (|6.4p and (|6.5p into the equation for /, after some computation we obtain 
the following equation (where we group the terms of the same order) 

= /f - fL + /°(P°)[V-a - W,] 

+ {gt - Lg) (6.6) 

- N{g), 

for 

N{g) =Xo{R^^Rt) 

+ Rl i^'pO-a - S'pO+a) 

- 9x(.P°) [{^p-a - V-a) - {Sp+a " V+a)] (6-7) 

— i?2 {Sp^a — Sp+a) 

^■.N\g) + N^g) + N^g) + N^g). 

Thus we have proved that 

Proposition 6.2. Fixed any admissible equilibrium f^,p^, then problem (|2.ip - (|2.2p is equivalent to 

' gt - Lg + N{g) 

< g{x,0)^gi (6.8) 
_g,(-AO) -ff,(B,0) = 0, 

for some g{x,t) that satisfies 

f + g{;t)GV, 

where L and N are given in (j2.7p . (|6.7p . respectively, and gj = fi — ■ The equivalence is simply 

f = f+9. 
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6.1 Lipschitz estimates 

We have that the non-Hnear part N is locally Lipschitz in Up-. 

Lemma 6.3. Fix p < c \^ for some c positive constant less than one. We have that the operator N 
given by ()6.7p is Lipschitz in the neighborhood Up, given in (|6.ip . More precisely, 

\\N{g)~N{g)\\x<C\\g^~g\\Y, 
for all g,g CzUp. The constant C depends on p and Aq. 

Proof: Let g,g £ Up, with roots e {p^ — v,p^ + v). Let p — p^ + q, p — p^ + q. From (j6.3|) we can 
estimate the difference I? — ?! by 





1 


< 






AO 




1 


< 






AO 




1 


< 






AO 



\q'q\- sup + + 



sup 



iW + i) - ~9{p' + 0\] 



If 

it holds 



\q-q\ < 



1 



;\\9-9\\ 



C°(p°-v,p°+v)- 



(1 -c)AO' 

Now we compute \\N{g) — N{g)\\-^. The definition of negative fractional Sobolev norm of order —0, 
61 > 0, is 

ll^llff-e = sup 

In our case, we will take 9 — — r + 2q: > 0. Note that for an exponent 3/2 < 6* < 2, we have that 
H'^{-A,B)cC^'^{-A,B). 

Since = bp±a — 5pO±a — Q^'pO^a^ '^^ easily estimate ||7V^(g) — A^^(5)||^, as follows: 



- ^ < ll'^pia - Sp±a\\H-0(^-A,B) + \l ^ ' PpO h-<> {-A,B) ' 

The first term can be bound by 

\\Sp±a - Sp±a\\H-e(_A B) - ^up { ± a) - ± a) | } < |p - p| = |g - g| , 

taking into account that H^^^ C L°° for 3/2 < 6* < 2. 
For the estimate of N^{g) - N^{g), it holds 



< 



X 



where 



< 1.9b) - ~9ip)\ + \9ip) - m\ + W) - 9ip")\ 

< 2 \\g — .9|lco(pO_^ pO+t,) + \\9\\c^p«-v,p'>+u) 
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On the other hand, for N^, 

\\9x{p") [Sp-a - SpO^a] ~ 9x{p") [Sp-a " 6pO~a] \\ x 

< \9iP°) -9{P°)\ Pp-a - V-allff-e + \9x{p°)\ \\Sp-a " Sp-aWn-e 

< 11.9 ^ .9llcO(pO-i/,pO + iy) + \\9\\c^pO-u,pO+u) W^p-a ^ <5p-a||^-8 • 

And for the last term N'^, 



RlSp-a — R2Sp-a < 11^2 {Sp-a — Sp-a) 



X 



^i?2 — -^2^ Sp-a 
R2 — R2 



Since 
it holds 



< |i?2| \\Sp-a - Sp-a\\H-o + 

I-R2I = \9x{p) - 9xip°)\ < \\9\\cHpO-„.pO+^) \q\ 



X 

\\Sp-a\\H-<> 



i?2 ^ R2 



< \9x{p) - 9x{P)\ + \9x{p) - gx{P)\ + \gx{p") - 9x{p°)\ 

^ Il5llc2(p0_^,p0+^) \q~ q\ + 2 \\g - 9\\c^p"-i^,p"+i^) 
< 



+ 2 \\g-~g\\y. 



,A"(l-c) 

In particular, there exists a constant C, such that 

\\Nig)~Nmx<C\\9~9\\Y, 
for every g, g e Up, where C = max |l, xo{i-c) } ■ 



(6.9) 



□ 



The previous estimate can be improved near the origin, in fact, DN{0) — 0. This estimate is needed 
in order to pass from local to global existence. More precisely. 



Lemma 6.4. We have that 



=o(ll.9lly), when \\g\\ 



0. 



(6.10) 



Proof: The previous lemma implies that ||-^(5)||x ^ O (WgWy) when WgWy ~^ 0. For the improvement 
we need, take a perturbation of the form f = + £g for some e > small. Let p be the root of / near 
Po. Because of (j6.3p we can write p = p^ + eq for q — Assume ||.g||y = 1, for simplicity. First, note 
that because of the Taylor theorem we can write 



Rf - 



^xt±a^^'p°±a 



for some 1 0:3 — | < eiy. We compute 

S'pO_a\\^.s ^ sup |^'(a;3-a)-(/)'(/-a) 



(6.11) 



where 9 = —r + 2a. Then, because of our choice of 9, H^{—A, B) C C^'^{~A, B), and from (|6.1ip we 
deduce that 



As a consequence. 



N\eg)\\^<Ce'+^. 
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Next, note that Ri — e [g{p) — g{p'^)] , so that |_Ri| < \\g\\ci ^-i^d this gives a bound for N"^. We 
handle the rest of the terms in a similar manner, and we conclude that 

\\N{g)\\^ < C,e'+^. 

for some C2 depending on (Yzt^att, but independent of The lemma is proved. □ 

Remark: The choice of s in the definition of the function space Y is made here. In particular, we need 
(j2,;3 regularity near the origin in order to have a Lipschitz estimate in Lemma 16.31 



6.2 Local existence 

First we review some basic facts about the nonlinear evolution equation 

i,.-Lg + Nig), 
\5(0)=5/- 

Here we always assume that L is a sectorial operator, U an open subset oi Y = X" for some < 
a < 1 and N : U ^ X & locally Lipschitz function. More precisely, given any g ^ U, there exists a 
neighborhood V C U oi g such that for all g,g ^ V, it holds 

\\Nig)^NCg)\\<C\\g-~g\\y, 

for some constant C depending on the neighborhood V. 

A solution of this initial value problem on (0,ti) is a continuous function g : [0,ti) X such that 
g(0) = gj- on {0,ti) we have g{t) e U, g{t) G D{L) = Z, ^ exists, t N{g{t)) is locally Holder 

continuous, /g°^^ \\^i9{t))\\ dt < 00 for some > 0; and moreover the differential equation (|6.12p is 
satisfied on (0, ti). In particular, we have the Duhamel's formula 

.Lt^. f\L{t-s). 



git) = e^'gi + e^('-^> N{g{s))ds. 



It is well known that 



Theorem 6.5 (Theorem 3.3.3 in 12J). Assume that L is a sectorial operator, < a < 1, and N : U ^ 

X , U an open subset of X°', and N is locally Lipschitz. Then for any gj € U, there exists < T = T[gi) 
such that (|6.12[) has a unique solution g on (0,T) with initial value g{Q) — gi. 

Remark: In fact, it is true that 

5GCi([0,T),X). 



Let gi = fi — g^ ■ The previous theorem assures local existence in time for g \i gj G Up. Then, 
Proposition lG. 21 gives that, the function f{x, t) = f^{x) +g{x, t) we construct is a solution to our original 
problem. Moreover, it has the right signs: it has a unique root in [—A, B], call it p{t), and it is positive 
in [—A,p{t)), negative in {p{t),B\. The proof is similar as in Lemma 
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6.3 Existence of a center manifold 

In the present section we give the proof the existence of a center manifold for our problem. Here we use 
follow the approach of chapter 6 in the book [12', that we find very clear and suitable for our purposes; 
although is the classical one for applications in partial differential equations. Other references are, 
for instance, [7] and [2]. 

We consider the equation (|6.12p . A set S C X" is a local invariant manifold if for any gj € S, 
there exists a solution g{-) of the differential equation on an open interval (ti,i2) containing with 
g{0) = gi and (/(•) e S for ti < t < t2. We say that S is an invariant manifold if we can always choose 
(ii,i2) = (-oo,+oo). 

In particular, for a coupled system 

i (5i)t ^1(51,. 92), 

\ (.92)4 = ^2(51 + 92) + N2{gi,g2), 

we look for invariant manifolds of the form S = {(.91, .92) '■ 92 = ^'(ffi)}- It is well known that: 

Theorem 6.6 (Theorem 6.1.2. in [12j). Let Xi,X2 he Banach spaces and assume —L2 is sectorial in 
X2. Let U be a neighborhood of the origin in X2 for some a < 1. Suppose that 

L The operator N2 '■ Xi x [/ — > X2 is locally Lipschitz with 

\\N2i9i.92) - N2{h,~92)\\x, < + \\92 - Mx^) ^ (6-13) 

and 

||A^2(gi,52)|lx<^- (6-14) 

2. The operator Ni : Xi x U Xi is locally Lipschitz with 

11^1(31,52) -A^i (31, ff2) II </^ll3i-3i|lxi +^'^2||32-.92|lx,- (6-15) 

3. There exists some 70 > such that 

||e^^*|L < Cie-''"*, ||(-L2)"e-^=*|L < Cii-"e-T°*, for t > 0. 

4- For some positive constant D, 

[92 ■■ \\92\\x^ < C L/ and CiN t-^e'^^^dt < D. 
5. For some positive constant E, 

poo 

e := kCi / w-"e-^°"e('^+^*'^^)"dw (6.16) 



has 

E f (1 + E)M2 1 

e< and 61 max U,-^ ^ ^ < 1. 

- l + E I M + M2E J 

Then there exists an invariant manifold (the center manifold) 

S = {(51,52) : 52 = cr(5i), -00 < t < +00,51 e Xi} , 

with 

\H9i)\\xs<D, 

and 

lk(5i) - ct(5i)IIv° ^ ^ 11.91 - .9i|lxi • 
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Remark: If /i < 70 and k is sufficiently small, we can easily satisfy condition 5 above. 
Remark: Note that condition 4 a-nd 5 ensure global in time existence of solutions, that have been 
proven to exist only locally in time in Theorem 16.51 In fact, if the constants Ci, C2,, N etc satisfy the 
above conditions, it is easy to check that (see Theorem 5.1.1 in 

ll.92|lx=<P, Vte(0,T). 

This implies that the solution exists for all time (Theorem 3.3.4 in [T^]) and stays in the neighborhood 
U, i.e., 

\\92\\x?<P^ Vte(0,+oo). 



In the following, we check that Proposition 16.61 can be applied in our problem. First, Proposition 
15.51 shows that we can split X — Xi © X2 with Xi = ker(L). The subspace Xi has finite dimension 
two. We write g ^ gi + 92, gi G -^1 and g2 G X2- Let Li be the restriction of the linear operator L to 
Xi, i = 1,2, and let P be the projection on Xi along X2. We denote 

Ni{gi,g2) ■=PN{gi+g2), N2{gi,g2) := [I - P)Nigi + g2). 
Thus, our problem (|6.12p is equivalent to 

({gi)t=Ni{gi,g2), 

\ (52)4 = -^2(32) +^^2(51,. 92)- 

In order to apply Theorem l6.61 we need to check the Lipschitz condition of the non- linear terms. In the 
following proposition we give a characterization of A^i and A^2- 

Proposition 6.7. Consider any g{-,t) G X such that /° + g{-,t) has a unique root in [—A,B], denote 
it by p{t). Then there exist constants 71 and ^2 depending only on A,B,a such that if g is written as 
g — gi + 52, gi G Xi, i = 1,2, we can write 

Ni{gi,g2) = 7i^2(i)go + l2R2{t)ho, 

where gQ, hg are the eigenf unctions of the zero eigenvalue given in Proposition \4-^ and R2{t) is defined 
as in (|6.4p 

R2{t) -.^ gx{p{t)) - gx{po)- 

Proof: We write 

Nig)^N,{g,,g2) + N2{gi,g2), (6.18) 
where Ni{gi,g2) G Xi, i = 1,2. Because {g^, /iq} is a basis for Xi, we can write 

Ni{gi,g2)^k{t)go + z{t)ha, / ^^2(51, 52) da; = / N2{gi, g2) dx = 0, (6.19) 

J -A Jo 

for some k{t), z{t) functions. The exact value of these coefficients can be computed thanks to Lemma 

EH 

On the other hand, consider the expression for N(g) given in (|6.7p . If we integrate it from to B, 
for each fixed time, we obtain that most of the terms cancel and we are left with 

hWig)] := [ N{g)dx^R2{t). 
Jo 

A similar integration from —A to gives that 

Ii[N{g)] := r N{g)dx = -R2{t), 

J -A 
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where R2{t) is defined as in It follows that Ii[N{g)] = -l2[N{g)] and from ([O)) we get 

2(-a + A + B) A-B 
a(a — 2A)[a — 2B) (a — 2A)(a — 2B) 

as desired. □ 

With all these ingredients we can prove that: 

Proposition 6.8. Fixed an equilibrium point , there exists some p > such that there exists a center 
manifold for (j6.17p in the neighborhood Up with 

lk(5i)IL<A 

and Lipschitz condition 
for D, E given as in Theorem\6.i 



Proof: It is a consequence of Theorem 16. 6[ but we need to check carefully the constants involved. 
First, note that the operator —L2 generates an analytic semigroup and satisfies the estimates given in 
(1521) 

l|e^'*llv <Cie-'"'*, ||(-i2)"e-^^*|L < Cii^^e-'^''*, 

for some < a < 1 and < 70 < 7. 

Next, we look at the non-linear part. We will work in a neighborhood Up for some p < min{Ao, vXa} 
so that we are in the situation of Lemma 16.11 

Lipschitz condition (|6.13p follows from Lemma [^31 On the other hand, (|6.15p is a direct consequence 
of Proposition 16.71 and (|6.9p . The Lipschitz constants K,p,Al2 depend on p. More precisely, these 
constants can be estimated by 

J^^Oip) in Up, 

when p < cA° for some c < 1. Fix p^ such that condition 5 is satisfied for all p < p5. 

Now, in order to check condition 4, we need Lemma I6.4[ since the Lipschitz estimate for N is not 
enough. We have seen that in this case, 

\\N{g)\\^^C2\\g\\'y+^ , when||5||y^0, 
for some < /3 < 1 and C2 = O ( %ii ) • Then looking at (|6.14p we can write 



(l-c)A" 

2 



iV=^pi+'3 in Up. 



Choose D — cp , with c < 1, we get 



Choosing p small enough, condition 4. is satisfied. 

We conclude that the solution exists globally in time and there exists an invariant manifold 

S = {(51,52) : 52 = cr(5i), -00 < i < +00,51 e Xi} , 
as we wished. □ 
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6.4 Stability of a family of equilibria 

Now we try to understand the decay near the center manifold. We are given an initial condition // 
with masses 

rPi fB 

fi = mi, - / // = 1712- 

-A Jpi 

Assume that it satisfies ||// ^ /''Hy < P for some p small enough that will be specified later. Here 
is any of the stationary states described in Section [3] with masses mi , m2 (not necessarily equal to 

mi, 1712)- We can take, without loss of generality, p° = 0. 

By linearizing around f'^, we are going to show that there exists a solution / of the system with 

initial condition //, that exists for all time, and that decays exponentially to a stationary state foo- 

Moreover, foo is chosen as the unique stationary state with 

Boo l-B 

foo = mi, - / /oo = m2, 

A Jp^ 

for /ooboo) = 0. 

Global existence and asymptotic decay are given by 

Theorem 6.9 (Theorem 6.1.4 and Corollary 6.1.5. in [12'). In addition to the hypothesis of Theorem 
\6.6l assume that dimXi < 00, and 

r:=e(l + M2 ^^^} \< 1, (6.20) 
V 70 - M / 

for /i' — 11 + M2E. Then S is uniformly asymptotically stable with asymptotic phase. Specifically, there 
exist S > 0, c > so that any solution {gi{t) , g2{t)) with initial condition satisfying 

11.92(0) -a(gi(0))||^. <<5, 

exists for all time t > and there is a solution gi (t) of 

{gi)t = Ni{gi,a{gi)) 

such that, for t > 0, 

hiit) - 51 (Oil + \\92{t) - <7(-9i{tmxs < Ce-'^'^ IMO) - <j{9im\\j,. , 

where ^ 

7 = 70 - (7 - n')r~ > 0. 

Remark: In the proof of this theorem, 6 is chosen so that the solution stays in the domain of existence 
for all time t > 0, that it is equivalent to 

kCie-^^*'>S < D/2. (6.21) 

Moreover, 

^ ^ ;~' 

where k and K are a constants depending only on a. 
Thus in our problem. 

Corollary 6.10. The center manifold constructed in Provosition 1 6. 8\ is an attractor if we start with 
initial condition gj G Up for some p small enough. 
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Proof: We take p as in the proof of Proposition 16.81 The value of 9 is given in (|6.16p . It is clear that 
we can take a smaller p so that (|6.20p is satisfied and p < 6 ioi 6 given in (|6.2ip . The corollary is a 
consequence of Theorem 16.91 and the remark afterwards. □ 



Up to this point, we have reduced the problem to understand how the trajectories look like inside 
the center manifold. 



Proposition 6.11. Any trajectory in the center manifold of Provosition wTS] must be stationary. 
Proof: A trajectory in the center manifold is given by 

5 = 51+52, gi = c{t)gQ + d{t)hQ, g2 = cr{gi). 
In order to find the coefficients c{t), d{t), we substitute this solution in (|6.17p . 

ego + dho + dt {a{gi)) = L{gi + gs) + ^(51 + 52), 
and we project onto Xi, thus 

ego + dho = Ni{gi +52)- 

We have an explicit formula for the right hand side, thanks to Proposition 16.71 Then, looking at each 
eigenspace, we obtain a system of ordinary differential equations for c{t), d{t): 



(c = JiR2{t), 

\c? = 72i?2(i), 

where 71, 72 are explicit constants depending only on a. A, B. 

We claim that this system of ODE's has a unique global solution. For this, consider 

t) = c{t)go + d(t)ho + <7{c{t)go + d(t)ho). 

It is easy to check that, for —v< p{t) < v and go, ho as in Proposition l4.21 it holds 

gM),t) = + f^U^), 5.(0, t) = c(0(l + 

Taking into account that i?2(t) = gx{p{t),t) — gx{0,t), the above system becomes 

\d = 72dW(af_^,,-<.^J. 



(6.22) 



(6.23) 



Theorem 16.61 ensures that a is a Lipschitz function. It follows that the right-hand side of (|6.23p is also 
a Lipschitz function with respect to the function c{t) . Existence of a unique solution locally in time is 
proven. 

We claim now that c{t) — cj and d{t) = dj solve system (|6.23p . Let gi{x,t) — cigo{x) + diho{x). 
We integrate for x £ {—A, 0) equation 

(5i)t = ^^1(51,52), 

we obtain that 

d f° f° 
0=—/ gidx^ Nidx^R2{t). 
dtJ^A' J -A 

Thus, the claim is proved. This implies that system (|6.22p has an unique global in time solution 
c{t) — ci, d{t) — di (the constant one). □ 
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In general, the center manifold is not unique and it does not consist only of the stationary points. 
However, in our case, the set of stationary states has dimension two and is contained in the center 
manifold. This is a special case and appears in other works, see for instance [9j on the Navier-Stokes 
equation, that has a similar behavior. In particular. 

Corollary 6.12. Our center manifold is exactly the set of equilibrium points for (|6.17p . i.e. 

{geUp : g = 91+92,91 Xi,g2 = cf{gi)} = {g e Up : Lg + N{g) = 0}. 



Now we can finish the proof of Theorem 12.21 For each initial datum // in the neighborhood Up, 
we have built a solution of the nonlinear problem for all time t > 0, of the form f = f'^ + g, and that 
converges exponentially to an stationary state, that may not be the same as f^ . Since we have already 
classified all the stationary solutions, in particular, because of conservation of mass, it has to be the 
precise /oo- 

Note that any stationary solution /oo can be written as a shifted zero eigenvector in the following 
way: 

>^oo9o{x ~ Poo) X e [-A + poo,B], 
AooSol-^) X e[-A,-A + poo]. 



foo{x) 



6.5 Some final remarks 

In this last section, we group together a few remarks. First, we can write a more explicit formula for 
the center manifold. Consider /° a general equilibrium point with /'^(O) = 0. Consider /oo the unique 
stationary state corresponding to the initial datum f^ + gi, with gj eU p. Let gj be decomposed into 
our basis of eigenfunctions 

gi = ciga + diha+ gi, gi e X2. 

The corresponding solution f{x, t) = /"^ + 51(2;, t) + g2{x, t) converges exponentially fast to the equilib- 
rium point 

/oo = /° + c/go + diho + a-{ci9a + djha), 
where cjgo + diho + cr{ci9o + djho) belongs to the center manifold. It holds that 

hi{ci, di) := / foodx = Aoofl {p^ - ^ + A) - = (c/ + A°)/i[go] + d//i[/io], 

J -A 

h2{ci,di):= / /ooda; = - Aooa(-B-poo - f) = (c/ + A°)/2M +d/-f2[M: 

Jo 

taking into account that any general equilibrium point /° can be written as f'^ — X'^go- Inverse function 
theorem shows that given c/ and d/, one can uniquely determine /oo (or, equivalently, determine poo 
and Aoo) from the above formula. More precisely 

(Aoo,Poo) = H^^{hi[ci,di),h2{ci,di)), 

where 



_ X„2 

H{\,p) 



Aa (p - f + A) - \p 
^p^^Xa{B-p-^ 

In fact, we can compute its derivative DH{X'^,0) and get 



DH{X'',0) 

that is positive definite because 



-af-f+B) X°a 



det DH{X°,0) = X"a^{B + A - a) > 0. 



29 



The inverse function theorem states the existence of a neighborhood near (A°, 0) such that the function 
H is invertible (and continuous). Therefore one can write our function a as 

CT(c/.go + diho) ^ -(c/ + A°).go + c?/^o + foo{H^^{hi{ci,di), h2{ci,dj))). 

Remark: The regularity of the free boundary p(t) is given by the following argument: since f{-,t) e 
C'^'^{—v, u) and fx{x,t) > cA° for all x G (— i^, i^), the function 

,f,. fxx{p{t)) 

' = -urn 

is uniformly bounded in time. This implies p{t) is Lipschitz regular \/t > 0. 



And to finish, we would like to see that the size of the neighborhood does not depend on the 
equilibrium we start with. This is achieved through a careful control of constants. In fact taking as 
family of equilibria all the functions such that 



A := {/"admissible : A° > x} , 



we can see that the constants C and C2 in Lemma 16.31 and 16.41 are uniformly bounded with respect to 
A*^. Tracing back the constants in the proofs of this section, we have Theorem [ 



7 Appendix 

In this last section we give the proof of Lemma 12.11 The proof of this result in a more general setting 
is included in the forthcoming paper For completeness, we briefly summarize here the main steps. 
Assume, without loss of generality, that A = B = 1, a = 3/4. 

Consider the Hilbert space HP(T^) for any p G M, with norm given by 



(l + n2)f|/„|2, 



where /„ are the Fourier coefficients of / in [—1,1], i.e., f ^ J2 fnen- 

The second functional space we deal with is H^{T^), which is defined through the norm 

wnH: ■■=uf\\m- 

Here the function (j) denotes a smooth cutoff function with support on [—1/2,1/2], identically one on 
[—1/4, 1/4], strictly increasing on [—1/2, —1/4], and strictly decreasing on [1/4, 1/2]. The intersection 
space X := n iJ' has a norm defined as 

ll/llx ■■=\\f\\Hp + \W\\m- 
The Laplacian operator acting on X can be written as, 

A/ = ^ A„/„e„, A„ ^ -n^, 

where A„ are the eigenvalues of the Laplacian. 

Our aim here is to find a functional space Z dense in X such that A : Z ^ X is a bounded operator. 
We claim that this functional space is given by 





30 



where (j) is another cutoff function such that = 1 on the support of 0. 
The choice of the cutoff functions <j) and cf) implies that 

0A(0/) = (/)A(/). 

As a consequence, the operator A : Z ^ X is bounded. 

The proof of the density oi Z in X requires some technicahties: first, we will prove that periodic 
C°°-functions are dense in X . This means that given fo € X and e > 0, we will find a smooth sequence 
fe such that II /e ~ faWx < 

By taking a collection C of overlapping intervals covering [—1, 1], we construct a good approximation 
for /o in each subintervals of C Consider for example 

C = {[-1, -1/2], [-5/8, -1/4], [-3/8, 3/8], [1/4, 5/8], [-1/2, 1]}. 

Once convergence in each single subinterval is shown, convergence in the whole interval [—1, 1] follows 
from the convergence in the single subintervals by a partition of unity argument. 

The convergence in is guaranteed by the regularity of the heat equation. Indeed, let f{x,t) be the 
solution of the heat equation ft = Af with initial condition /q. Then, if /o = fn^m it is well known 
that 

satisfies / ^ /o in when t ^ 0. Moreover, also / ^ /o in i?^ on [—3/8,3/8], away from the 
vanishing points of the cutoff 0. 

On the subinterval [—5/8, —1/4] we regularize the function /o as fe — d^* /q where is a mollifier with 
support in (0, 1) and 9^{x) — -6{x/e), with supp 6^ — (0, e). Since /o G W, it is known that /o 
in HP. 

It remains to prove that (t>fe converges to 0/o in H'^: since (/)/o G H'^^ one can consider a small 
shifting (t){x — t)fo{x) E H'^ for t £ (0,e), so that <j){x)fQ{x + t) E H'^. In this way, the function 
4>{x)fe{x), defined as 

=0(x) / U{x + t)e,{t)dt, (7.24) 



still belongs to iJ* with norm uniformly bounded with respect to e, depending on H^/oH/f, ■ 
The term \\(j)ft — <t>fo\\Hi '^^'^ rewritten as 

110/. - <i>h\\m < Uife - A)\\h^ + MA - n\H^ + Uih - iM„, , 

where ^p^ := 9e * ip- 

The first term of the sum can be uniformly bounded by \\(j){fo — < '^/^ (in the same way (I7.24p 

is bounded by ||0/o||j:^,). 

Due to a standard property of mollifiers, we have that ||0(V'e — '>P)\\hi — ^f^- 

Finally, for all 6 > 0, there exists a function ip G C°° with compact support in the interval [—5/8, —1/4] 
suchthat ||0(/o-V)llff, <'5/3. 

This implies that for every e small enough, one can find a S such that ||0/e — ^/o 11^-9 < S. 

Once the proposition is true for the A operator, then it is true in our case, where we deal with a 
bounded perturbation of the Laplacian: 



Lg -.^ gxx - gx{0) 6^^_3 - S^^;> +5(0) 
It is easy to check that 



S' 3 - S' 



\\Lg\\x < l|A5llx + ll.9llcl <^ll5llz- 

□ 
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